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1. INTRODUCTION 



1,1 The purpose of this,. program is to produce data for gain-phase, Bode and 
Nyquist plots from transfer functions in the s-plane. The operator may- 
supply the computer with either un -normalized polynomials or normalized 
factors, 

2. METHOD 

2.1 The program produces data for frequency response plots from transfer func- 
tions having the form of: 
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2.2 If the transfer function is entered in factor form, it is first expanded 
into a polynomial. Then, for entered values of w the transfer function 
G(s) is evaluated to determine the imaginary portion of G(jw), the real 
portion of G(jw), |G(jw)| , 20 log 1Q | G(jw)| and the phase angle of G(jw). 
These results are then printed out. 

3. RESTRICTIONS 

The maximum order of the numerator and denominator s-plane polynomials is £6. 

The maximum number of linear s-plane factors is £6 and the maximum number 
of complex s-plane factor pairs is 28. The total number of s-plane factors 
(linear and complex) should not exceed S>6. 

k. USAGE 

U.1.0 After loading the tape, put sense switch C down if the form of the transfer 
function is normalized factors, and proceed with step U«2«l-' If the form of 
the transfer function is polynomials, put sense switch C up and proceed with 
step li.l.l. 
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U.l.l Press the Start 1 button on the right side of the control console. The 
typewriter will type an N. 

U.1.2 Enter the degree of the numerator, NN, by typing two digits from "00" to 
"60" • The typewriter will then type A(NN) where NN is the degree of the 
numerator just typed. 

U.l»3 Enter A- in floating point from the typewriter (see floating point input 

instructions lu3«0. The typewriter will then type A(NN-l) after which A- T , 
must be entered, etc. until A(00) is typed and A^ is entered. 

U.l.Ij The typewriter will type M. Enter the order of the denominator in the same 
way as the order of the numerator was entered, and enter the denominator 
polynomial in the same way the numerator polynomial was entered (See ii.1.2 
and U.1.3). 

U.1.5 The typewriter will type R. Enter the number of zeros minus pole s at th e 

origin by typing a sign (+ or -) and two digits (00 through 99) • ^ ~^ 

li.l«6 The typewriter will type K. Enter in floating point the transfer function 
Main. (See floating point input instructions, Section U.3.0) . 

4.1.7 - The typewriter will type "ERROR?". If an error was made during any of the 

above input procedure, type those characters typed by the computer preceding 
the error and then the correct type in. If no errors were made, or after all 
errors have been corrected, type a carriage return. 

I4.I.8 The typewriter will type W. Enter in floating point the frequency. in radians 
at which the transfer function is to be evaluated. (See Section U.3.Q) • 

^.1.9 The typewriter will type DW. Enter in floating point the incremental frequency 
in radians which you wish W to be incremented by for each evaluation. 

lwl.10 If a new W and DW is desired, set sense switch B down and W will again print 
out. Set sense switch B up and proceed from ij.1.8. 

U.2.1 Press the Start 1 button on the right side of the control console. ,The type- 
writer will type C. Enter the number of linear factors in the numerator by 
typing two digits ranging from 00 to 60. The typewriter will type E(CC) where 
CC is the number of linear factors in the numerator. 

ii.2.2 Type the value of E c in floating point. (See Section U.3.0). The typewriter 
will type E(CC-l) . Type E c , and continue this procedure until all the linear 
numerator factors have been""entered. The typewriter will then type D. 

U.2.3 Enter the number of complex pair s in the numerator by typing two digits, DD. 
The typewriter will then type F(DD). 

U.2.U Enter in floating point the real part of the first complex pair (See Section 
ii.3.0). The typewriter will then type G(DD). Enter in floating point the 
imaginary part of the first complex pair. The typewriter will then type 
F(DD-l) etc. until all of the numerator complex factors are entered. 
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ii,2.^ Enter the denominator in the same way as the numerator entering the number 
of linear factors after P is typed, the linear factors in floating point 
following H(PP), H(PP-l), . . . HOI, the number of quadratic pairs after Q 
is typed, the real part of the quadratic pairs after l(QQ), I(QQ-l), . . . , 
101, and the imaginary part of the quadratic pairs after J(QQ), j(QQ-l), 
. . . , J01. Then follow statements U.l.£ through U.1.10. 

U.3.0 Floating Point Input 

In entering floating point numbers, tyoe the sign of the mantissa, the man- 
tissa assuming the decimal point is in front of the first character typed 
(DO NOT TYPE THE DECIMAL POINT), the sign of the exponent, the exponent, and 
a carriage return. Examples: 

Number Type 

1 +1+1 (carriage return) 

-10.76132 -1076132+2 (carriage return) 

.75 x 10" +7^-26 (carriage return) 

+0+0 (carriage return) 

2 U0 _ 2 U0 

Numbers may range between 2 and 2 and any number of digits may be 
typed for the mantissa, though the computer* will accept no more than approxi- 
mately 10 place accuracy. 

U.U.I If a transfer function with a polynomial numerator and factored denominator 

is to be entered, change the position of sense switch C af ter(W has typed out. 

U.U.2 If a transfer function with a factor numerator and a polynomial denominator 

is to be entered, change the position of sense switch C after C has been typed. 

h.5.0 Program Results 

After computing the frequency response, the computer will type six floating 
point numbers which are the values of: the frequency in radians, the imagin- 
ary part of G(jw), the real part of G(jw) , J G( jw)| , 20 log 1Q | G(jw)| , and 
the phase angle in degrees. ' M 

U.6.0 If an error is detected at any time after Section U.1.7 is completed, stop 
the computer and press Start 2 button on the right of the control console 
and proceed with Section U.1.7. 

U.7.0 To. enter a new problem, stop the computer and proceed with Section U.1.0. 

5.0 EXAMPLES 

See Appendix. 
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6.0 CODING INFORMATION 
6,1,0 Storage Requirements 

Program Storage 0000-2120 

Working Storage 3000-3170 

6.20 Timing 

Computation time varies with the size of polynomial. The time between 
successive frequency response print-outs is approximately 3 seconds. 



5.0 Examples 



to 



EXAMPLE NO. 1 



Transfer Function 216S' 



S 6 + 9S 5 + 



S^T 



21S^ + 



3QS k + 



90S 3 + 12!iS 2 + 96S 



32 



208S' 



1190S" 



UlOi^S^ + 806US + 6912 



n o6ao6 +i+ i 
A05 +9+1 

AW +38+2 
AO3 +9+2 
A 02 +12H+3 
A01 +96+2 
A00 +32+2 

M O6BO6 + 1 + 1 
BOS +21+2 
804 +208+3 
BO3 +II9+U 
B02 +k)Ch+k 
301 +8054+^ 
BOO +69I 2+H 

R +02 
K +21 6+? 
ERROR? 

W +1 + 1 

nw +1+1 



Sense Switch "C M is up. 



W 



Frequency 

09999999 00 
19999999 01 
29999999 01 



Imag* G(jw) 

-161891^3 01 
-86cn 90*19 01 

-28578^61 02 



Real G(jw) 



G(jw) 



201g__ G(jw) 



3 10 







23800286 00 I6363166 01 U2773U72 01 -8163629I 02 
21111312 02 22796119 02 27157218 02 -22166251 02 
I7IO2769 03 17339896 03 IA780929 02 -91*863505 01 



EXAMPLE NO. 2 



C 02 
E02 +1+1 
EOI +2+1 

D 02 
F02 +1 + 1 
G02 +1 + 1 
F01 +2+1 
G01 +2+1 

P 02 
H02 +3+1 
HOI +4+1 

Q 02 
102 +3+1 
J02 +3+1 
101 +4+1 
J01 +4+1 

R +02 

K +216+3 

ERROR? 
W +]+) 
DW +1 + 1 

Frequency 

99999999 oo 

19999999 01 

29999999 01 
w 



Transfer Function 216S' 



Imag. G(jw) 

-16189153 01 
-86009*9 01 

-28578^61 02 



(S + l) (S + 2) .(S + 1 A. jl) (S + 2 ±,J2) 

(S + 3) (S + 4) (S + 3 + 33) (S + 4 + 34) 



Sense Switch "C" is down. 



Real G(3w) 

23800286 00 

21111312 02 

I7IO276Q 03 



|G(jw)J 

16363166 01 

22796119 02 

17339896 03 



201og 10 |G(3w) 







^2773^72 01 -8163629I 02 
27157218 02 -22166251 02 
4^780929 02 -9^863505 01 



